


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 


DSpace Repository 


Theses and Dissertations 


1960 


1. Thesis and Dissertation Collection, all items 


Random neutron flux variations in a 
Subcritical assembly. 


Rough, Jimmie Lynn 


lowa State University 


Downloaded from NPS Archive: 


ai | pune 
mt \ LIBRARY 


http://www.nps.edu/library 





Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 

research materials and institutional publications created by the NPS community. 

Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholarly author. 


Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 


NPS ARCHIVE 


1960 
ROUGH, J. 





RANDOM NEUTRON FLUX VARIATIONS 
RP serine YT Cal 


JIMMIE LYNN ROUGH 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOO! 
MONTEREY CA vole 8-510 











Ree 








RANDOM NEUTRON FLUX VARTATIONS 
IN A SUSCRITICAL ASSIMBLY 


by 


Jimmie Lynn —_ 
] 


A Thesis Gupmitted tc the 
Graduate Faculty in Partial Fulfillment of 
fhe Requirements for the Degree of 


HASTIRN OF SCIENCE 


“ajor Subject: Nuclear Engineering 





DUDLEY KNOX | ARY 
Rei lAL POSTAR AH iA fe or 
a R wom 4.4 FE OF 
44 MONTEREY Cx dol @ G3 


TABLE OF CONTTATS 


Pare 

INTRODUCTION 1 
NEVIiW OF THE LITERATULE 3 
DEVELOPHERT OF THE KINETIC EQUATIONS 6 
DERIVATION OF THE TRANSFER PUNCTIONS 14 
Transfer Function with Source Variation 1h 
Transfer Function with Reactivity Variation 18 
Modulus of the Transfer Function a1 
STOCHASTIC EVALUATION OF THY TRANSFER FUNCTION 23 


DESCRIPTION OF APPARATUS AND EXPERIMINTAL PROCEDURE 29 


DISCUSSION OF RESULTS “1 
CONCLUSIONS LG 
SUGGESTICNS FOR FURTHER STUDY LF 
LITERATURE CITED LS 
ACKNOWLFDGHMTNTS 50 


APPENDIX STZ 





“ie 
“4 ome =) ae eed 
— eee caw ui 





@ettesas 95.577 7% 


INTRODUCT ION 


The subcritical assembly has been usec as an educational 
device to demonstrate the principles of nuclear selence and 
engineering, and for research in experimental reactor 
engineering. There are several advantages in the use of a 
subcritical assembiy. The assembly is small in size and 
relatively inexpensive in materials, construction, and re- 
quired apparatus. Cnly neminal shielding is required due to 
the low activity involved. Lattice parameters can be varied 
more easily than in a eritical assembly. 

The kinetic equations of the assembly are time dependent, 
non-linear differential equations which may be linearized in 
the derivation of the transfer function. The method of ex- 
perimental evaluation of the transfer function is to vary cne 
or more of the parameters in such a way that the limits within 
which the linearization is valid are not exceeded. This 
variation of the parameters is called “fercing". The usual 
methods cf foreing in a subcritical assembly are varying the 
strength of the neutron source or the reactivity in a periodic 
manner. The minimum periodic input signal to which the re- 
sponse can be measured is determined by the amplitude of the 
random or stochastic variations in the steady state flux 
Level due te the inherent statistical nature of the fission 
process. It is possible, in a system with large stochastic 


variations in the neutron flux, te exceed the range of line- 
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arity of the kinetie equations by making the magnitude of the 
input signal large enough to obtain a readable response. A 
technique has heen develepec, from studies of stochastic 
processes such as Brownian Motion, which makes it possible 

to evaluate the transfer function of an assembly. The 
stochastic or random variation in the flux is considered tc 
be the response of the system te stochastic scurce foreing 
due to the statistical nature of the basic phenomena invelved. 
This methed cf evaluation of the transfer function sheuld be 


least likely to invalidate the linear approximations made. 
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REVIEW OF THY LITERATURE 


A great amount of effort has been expended in the deriva- 
tion and evaluation of critical assembly transfer functions. 
One of the primary reasens was the need for an autematic 
control system to maintain a given power level in a reactor. 
The similerity ef the transfer funetion of a reactor to that 
of a servomechanism control system hed been noted, and it was 
determined that this type of system was suitable for reactor 
ecntrol. 

Owens, Crever, and Pigott (13) in 1949 proposed the con- 
cept cf treating a critical assembly as a "black box" to 
which a feedback control system eculd be applied in order to 
maintain a desired nower level. The proposal was based on 
the observation that the rocts of the kinetic equation of the 
assembly were negative. This indicated that the critical 
assembly transfer function was the same as the transfer func- 
tion fer a minimum phase shift network, and that the transfer 
eharacteristic could be represented by a Bode diagram. In 
that same year Franz (5) presented a more detailed derivation 
of the transfer function and showed that the standard tech- 
niques of servomechanism theory applied. Harrer, Boyar, and 
Krueoff (9) evaluated the transfer funetion of the CP-2 re- 
actor by sinusoidally varying the reactivity, and showed that 
the transfer function could be used in the design of servo 


loecps or regulating systems which include the reactor. 
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Meekup assemblies have been utilized fer the determina- 
tion of the transfer function cf proposed reactors. YUoland, 
Smith, and hice (3) measured the zero power transfer function 
of a fast critical assembly, ZrR-III mockup assembly of 
EBR~I, Mark III, by determining the change in power level due 
to oscillation cf the reactivity. 

HKefinement in the techniques of testing and the evalua- 
tion of transfer functions are ecnstantly sought, and the 
suberitical assembly provides a relatively simple vehicle 
for such investigations. A suberitical assembly was utilized 
by Axtmann, Dessauer, and Parkinson (1) for reactivity meas- 
urements and testing of pile compenents. The transfer func- 
tion of the Iowa Stote University subcritical assembly was 
initially evaluated by Fieci (14) utilizing the methed of 
source forcing. 

A new technique for evaluation ef the transfer function 
of a nuclear reacter cperating at a steady state by an analysis 
of the random, or stochastic, nature of the neutron population 
has been developed recently. Moore (11) has shewn mathemat- 
leally that the square of the modulus of the transfer function 
of a reactor is proportional te the Fourier transform of the 
autocorrelation function which is equal te the power spectral 
density function of the system. This was verified experi-+- 
mentally by Cohn (4) and by Griffin and Lundholm (8). Velez 


(16) evaluated the transfer function of the Ford Nuclear 
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Reactor and the subcritical assembly at the University of 
MNiehigan by the direct measurement of the autocorrelation 
funetion, but the results were reported te be inconclusive in 
the ease of the subcritical assembly. 411 references found 
describing the evaluation of the transfer function by measure- 
ment of the power spectral density function pertained to 


critical assemblies. 


—— i 


oe 


de 





- 
iin 


ne 


anal 








ec. 


DIVELCPYIND OF THD RINTTIc PQUATIONS 


The kinetic equations of the subcritical assembly can 
be derived from the general diffusion equation (6, p. 223 and 


74 ee LOL) 


avg. F p+5=Sb-i-2 (1) 


whieh is applicable to that region cf a subcritical assembly 
away from the extranecus neutron source where most of the 
neutrons are fission neutrens. The symbcls used are defined 
in the table of symbcls. The first term represents the icss 
of thermal neutrons due te leakage from the assembly; the 
second term is the less of thermal neutrons due to absorption; 
the third term is the neutron producticn or source terms and 
the right side of the equation is the non-equilibrium chance 
in neutron pepulation vith time. 

The source term is made up of prompt and delayed fission 
neutrons, and neutrons from an external scuree. It has been 
shown (7, p. 193) that the number of fast fissien neutrens 
produced is <2. a G, and further (6, p. 226) that | 
is the tctal fraction of delayeé neutrons the prompt neutron 


seurce term Sp is 


2T 


-B 
s, = (1 -/8 Ma 2.» Be (2) 


There are six knewn groups of delayed neutrons that 


account for 0.64 per cent cf the thermal neutren producticn. 
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these delayed neutrons accompany the radioactive decay of cer- 
tain fission products, called neutron precursors, which exist 
in an excited state. If & is the fraction of each delayed 
neutron group then the formation of the precursors of the 
group would be Ee? 9. When the concentration of a precursor 
group is denoted by Cy anc the decay constant as A,, the rate 
of decay would pe A,e,. The net rate of formation of pre- 
eursers to the i-th delayed neutron group would then be 
Cc K w 

oh * hi -_ 29 é X4e, (3) 
The total rate of preduction of delayed neutrons is > es 
and when multiplied by the resonance escape proba hi tty, RP; 
and the slowing down non leakage probability, ee the 


result is the celaved thermal neutron source term 


6 
o 
7 ki 
S, =pe c (4) 
a°P ce’ a2 oe 


Table 1. Table of symbols 


Be Buckling em. 

p Fraction of delayed neutrons 

B34 Fraction of i-th group of delayed neutrons 

C Number of delayed neutron precursors in the 
assembly 

Cy Number of the 1-th group of delayed neutron 


preeursers in the assembly 


eee Te yams CO Mamlan! Get Gee SommBees decals: candi 
IHkve Oke OTR VE, Mita! elie ,ersatety aolsal) vin 
ineethh ee HO te ee a PM dt ee 
ott Vb Oreererey ar Ve emtteenstt act me/ ~~ ore 
treme § to meetremee ee ee ASE Fee bie ne 
He 2 DoOtiertol Pr eee) Tee ext 'gosh Gt Aleey yeas = 


ao ped? Skee yor, serine Leyulel de-) at? of plowtae 


oe ak eh Fo 
pofigh 2 Of ccurthear Moda lal th mo eromtony Ve ney lerot aiff 
re] Sere Ce eee ee 
ant) ay fer EAU re Wabi teen aves yaballa wal ab 
~~ 9 qe eerts portine: Teeny bets tes we) at 


















»* ~ i» 
eu See a 
oa a ek ees | 
> o- : _-= uae 
ie of «6 aavttue teenie’ te mee 20NTR FQ 


— eon Geqale’ Yo ports AME Te poiseavt 5 
oO of enernmew Cordis Musial So ce 

. ad 
qernex vyetad be low Ais oe Oh tg} 
a 


=~ © =< Be — 





Lad 


fable 1. (Contirwe?} 





Cc 


“— 


o 3 


os Ww Ml SS 


ua Ss 


ree 











Coneentration of delayed neutron precursers 

Coneentration of i-th group of delayed 
neutron precursors 

shermal neutron diffusion ceefficient 

Frequency 

Effective multiplication factor 

Multiplication factor for an infinite assembly 

Reutron diffusion length in the moderater 

Prompt neutron lifetime in a finite assembly 

Prompt neutron lifetime in an infinite assembly 

Average neutron lifetime in 2 finite assembly 

Average decay constant of delayed neutron 
precursors 

Deeay censtant of the i-th sroup of delayed 
neutron precursors 

Total number of thermal neutrons in a finite 
assembly 

Number of thermal neutrons per unit velume 

Resonance escape probability 


Thermal neutren flux 


Tate of production of thermal neutrons per 
unit volume from all scurce 
total external thermal neutron seurce 


contribution to the assembly 
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Table 1. (Continwed) 
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~<L External thermal neutrcn source ccentribution 


neutrons 
per unit volume ain ie Bee 
ze Macrosccpie absorpticn cross section om. 7! 
c Time sec. 
2 
ae Fermi age cM. 
Vv Thermal neutron mean velocity em./sec. 
radians 
WW F — a 
requency ax. 
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The contribution cf the external thermal neutron scurce term 


will be assumed to be 
2 (5) 
Substitution of Equations 2, 4, and 5 into Equaticn 1 yields 
22. nv’. 5 9+ 0-2 )ke F Gee T 
qt a Te a 


4 (6) 
rpem 2 dye + tL 


® and c, are functions of space and time, and it can be shewn 
(6, p. 227) that the space and time variables are separable 
allowing the partial differential equation to be reduced to 
an ordinary differential equation. Consideration of the 


fundamental space mode only (6, p. 227 and 9, p. 33) yields 
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the equality 
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V6 = - BY (7) 
the diffusion length is defined by 
i = = (8) 
a 
Rearrangement of the equation yields 


2 
D=L2, (9) 


The mean life of thermal neutrens in an infinite medium is 


defined by 


Bee 


VW 2,¥ 


When rearranged this relation becomes 
2 = (11) 
a f 7 
o 


The flux, the effective multiplication factor in a finite 
medium, and the prompt neutron lifetime in a finite medium 


ere cefined by 


G = nv (12 
kk Bt 
oo eee (13) 
1 + L“B 
Me o _* (rh) 
1 + B 


The combination cf ikguations 13 and 14 yields 
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t, x 
-. = Be (15) 
@ 
The substitution of lquations 7, 9, 11, 12, and 15 into Iqua- 


tion 6 yields 


*» 


Qo. . (1°B? +1) B+ 1-6) 
. ja 2 A (16) 


5 & 
+p gw 2. rae + aL 
$7 


The combination of Equation 14 with Yauation 16 and rearrang- 
ing yields 

* 
Gn = [eC ~f)- 1] ae +p e7 2 Ayes+ “a (17) 
dt A * =. 
When Equations 11, 12, and 15 are substituted into «vation 


3 the result is 


de 
a= any he 14 
pe™ “ff 


Equations 17 and 18 are the ecupled kinetic equations ef the 
Suberitical assembly. 

It should be nected at this point that in arriving at 
Equation 17 the asswoption has been made thet the external 
source contributes a different number of thermal monoenerzetic 
neutrons to each unit volume of the subcritical assembly, 
the number decreasing with height in the assembly. The 
kinetic equations ean also be derived by considering the 


total neutron population rather than the neutren population 
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on a per unit volume basis. 

Consider that at a particular time the assembly contains 
ii neutrens, and the effeetive repreduction censtant, k, is 
less than one. The prompt neutron lifetime is ts the prempt 


neutron reproduction constant is k_, and the rate of decrease 


P 
of neutron population in the assembly is Cc, - 1) u/ p*, if 
neutrons are not supplied te the assembly from an external 
scurece. Absorption and leakage aceount for the decrease. 
If # is the fraction of delayed neutrons then the prompt 


neutron repreduction constant, k_, can be replaced by 


p 
(1 ~f). The rate of deerecase of neutron population with 


time is then civen by 


ais [ka - 2) - 1) w/* (19) 


cy 


Neutrons are supplied to the assembly by an external source 
and six croups of delayed neutrons. The total number of 
Gelayed neutrons of the i-th greup stered in the fissions 
products of the assembly at any time is Cy. if the decay 
constant of the nuclei containing the i-th sroup cf delayed 
neutrons is ras the neutrons of that group are emitted at 

a rate A,C, as fast neutrons. These fast delayed nevtrons 
must be corrected for fast non-leskage, ~~ >» and resenance 


escape probability, p. Summation of the terms in the non- 


equilibrium rate eguation fer the assembly ylelds 
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2 
= = Eze =f te 1] - + pe" a d,C4+ 8, (20) 


ihe total number of thermal neutrons being produced per 
unit time is a ef which a fraction 4 are delayed. the 
délayed neutrons are mot lost by leakage cr absorption until 
after they are emitted from the fission products, and hence 


the number of delayed neutrons stored in the fission preduects 





pe 
rium equation for the i-th greup of delayed neutrons is 


Le * , - 
is greater than Fs by the factor — 7 the none-equilib- 


ac yx i 


i. ee we 
ac = 3 a Ay Cy (21) 


pe 


Equations 20 and 21 are the same as Fauations 17 and 18 ex- 
cept that they are based upon the tetal number of neutrons 

in the assembly, and the external source term then accounts 
for all of the neutrons supplied to the assembly. The neutron 
flux at any point in the assembly can be considered to be 
preportional to the total nomber of neutrons in the assembly 
at a particular instant. The use of either set of ccurled 
kinetile equations in the derivation of the transfer function 


will yield the same result. 
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DERIVATICGN OF THEE TRANSFFE. FUNCTICHS 


the kinetic eguaticns and the resultant transfer function 
ean be varied by varying the scurce or reactivity in a pre- 
detervined manner. This variaticn is usually sinuscidal when 
purposely introduced. Nowever, the scurce strength and re- 
aetivity may vary even at "steady state" conditions. This 
is due te the fluctuation in the number ef neutrons available 
to the assembly caused by the natural statistical fluctua- 
tiens in the rates cf neutron abserpticn and fission. The 
resultant randcm variations in the neutren flux make possible 
the evalusaticn of the transfer function by proper analysis 


at steady state. 
Transfer Function with Source Variation 


For steady state conditions (using zero subscripts) 


Equations 17 and 18 become 


6 
: n 2 
G28 20 = [ki ~2)- 1] 3+ pe“) D 
- =0O= = @)- —z + pe a + (22) 
- ea 4 17 in ‘4, © 
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-B* T d , O 
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Consider a smsll variatien $< to be superimposed on -sZ@ 
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se that t= t+ S24 
nen, + Sn (24) 
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Substitution of Equation 24 inte Parnations 17 and 18 ylelds, 





(n, + eo 
dn . dSn . 8/0 oe - 
yet FD gcey, akahesies + § 2 (25) 
and 
dey aSe Cy; At (ny + Sn) 
“at = “at —s- WP es - A shox + Se,) (26) 
pe 


When the terms of Iquations 25 and 26 are expanded, the 
steady state comronents of Nquations 22 and 23 are collected 
and set equal te zero, the resultant equations are 


6 
as Ka) yy ot. -B°T 
Soe = Sn - Sin 4 pe 2. Ay Sey Se (27) 
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by transforming to Laplace notatien Dauaticons 27 and 26 become 
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sén(s) = ee Sn(s) - Pp Ses) + pe BT Dd gSeq(s) 


+ 3 (5) (29) 


sée,(s) = - A, Sn(s) = A, Se, (s) (30) 
pie 2° 


,earrangement of Equation 30 yields 


Se,(s) » ~, oey7 8 als) (31) 
5 fee T i 


,- 











(Peay © Ga 8 arate af4 A evi see Ge sutra one 


° et Pr. 2% m4 oft. 
Te =~ e o ' hae ® natal ae - 


-_ 
fae’ eve eatys = ast —: oN a 


od) ppt tengee ewe SS ink SY eeatrae’! fh pee an 
‘Meedte Liae ote fe tee ct ati aed Yo dept Wrarl ° 
ss pend Tmt tate foe 


vw sd egiys 25% ene SB - 02 tee 


ib 
(mr | | et: sae 


eee Mea, Saree POR eeblent of ow 
“Wptt 3 Oe ° foe seh - lale get. il 2 
ee sah a ' 


. C co. + — } ts) ah - 
tw oo — Med 7 isle — 


PLC OF ope sa 
ee 
ne 


= 






































16 


Substitution of Equation 31 in Equation 29 and collecting 


terms ylelds 


6 
rk 
JSn(s) po - Hh ~ wlan kre A - Sg) (32) 
L~ "AS oo Pe ’ 
By rearranging Iaquaticn 32, the transfer funetion becomes 
S via 
re (33) 
sas 6 {4 rs 
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—= 6 /3, pI 
By substitution of AS for > ~=—+- Equation 33 becones 
1 


a On) 
sf-k(1 -4) +1 -k Ate 


Multiplication of the numerator and denominator by (s +) 


and division by sf yields 


Sats} . yp (35) 
ss S$ + (A+ SA + As + Re - &) 


The denominater is of quadratic ferm, and when averace values 
of = 0.08, k = 0.5, tL = 0.001, and /* = 0.0064 are 
substituted in the coefficient of s, it is seen that X and 
ye are negligible compsred to i= so that 
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On(s) _ + 
6 ' 3° ” (LH) - d i a = 


The roots of the denominator are feund by use of the quadratic 


formula. 
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From the binomial expansion (a + b)*® =a + nb ..... where 
a 
b a, the radical expands to E - +=.) 2 ~ ave seen 
By substitution of the expandéd terms fer the radical, s 


= i 


becomes 

cttw fas [1-24] j 
First root 

gy = tg fr +1 - BE fe) 


second root 


22k 2iv b= j 
so = AH je 1 + RY mf > = et) + 





When A is considered negligible compared to the first term 
the resultant rect is 

b. & = (eee) 

2 t* 
the denominator is now seen to be factorable to a good ap- 


proximation into 


(9 +A Ms + R558) 
By substitution back into Equation 36 the result is 


2 2 dts —1 (37) 


1s = i 
5 (3 #X)(s + Fe) s + 44 


This is the final form of the transfer fumecticn. It shews 


that for a subcritical assembly one may ignore the effect of 
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the delayed neutrons in determining the change in neutron 
population due te a variation in the scvrrece. The terms con- 
taining the average delayed neutron cderay constant cancel 

eut, and thus the prompt neutron lifetime can be replaced by 
the average neutron lifetime. The assumption that the effect 
of the delayed neutrons can be ignored is valid for Koop 
0.995 (1). The subcritical assembly utilized in this in- 
vestigation has a K np of abcut 0.52, hence the above assurp- 
tion shovld be valid. ‘he same result, tquation 37, is ob- 
tained from iquations 17 and 18 when the delayed neutron terns 
econtaining r5e5 are deleted and the same derivaticn procedure 


as above followed. 
Transfer Punction with reactivity Variation 


Again utilizing Fauation 17 and Equation 18 for the 
transfer function and tquatiens 22 and 23 fer the steady 
state cendition, let a small variation ok be imposed upon 
k so that 

kK 


i 


k, + Jk 
n + Sn (38) 
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eC, =a, + Se 
i 1, i 
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where §k is the amplitude of the variation, and not the 
instantaneous value of the excess reactivity. Substitution 


of Equation 38 into Equations 17 and 16 yields 





ieee ee et 

ji oe om! .sereeee aff 42 Ge atOF © Oo ee oc leah 
imeem Rear Gane! et ee real et opens oe teed oad 
eee leet ot eee RANT Ad ret Poemed nt) ce Oe AOR 
oem eels SOM ei qnamte WAT comeertl: coriqe athe 
jy! Tul Aisew te brroeel a4 bee Goreires evudel cer Oe 

coh Ohad Gt MedSTERG citbewee Leet Bence ct Le) OPER 
ress Seed 4A) pend (EE-8 Tire In yas © eet one eee 
40 bh ARE antynnet yf lewne anes aet sbi ler oe Aunmnie pate 
memset acme Ueyzateh oaks eve Al bs anctoveny vey baal 
evethdern pOLepetyy) amas -c) how Mabecet Oem yiyh owoiwes 


. bomg ied arote oe 
— » @bhnetet witvtse® dele onli tee 


S-« a 
wei cacy Sf xabtnctt hue ¥2 sorrans! pelallbio elias” 


cried. hea 2 enegnyl PRS eas 
re fenel 00> maciitory Lines © eC eae aaa 














ah *,f*2 
mu pee ae 


-s i ov «eee , —T 


en yer Sa a eke a 
om pw bu 
“np eg ot <okee veuamtoanal 


me Se meena sinned owen 


«oc 4 == —_—" « 

















i> 





gn dons [ue + Sea ~~) - y | aes 
pe 
a A, a ~ Sa r* < (39) 
de ase (k= S$k)(n. + Sn) 
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By multiplication, setting preducts of small numbers equal 
te sero, and collecting the steady state terms of Equaticns 


22 and 23 and equating them to zero one obtains 
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ade, = Ak, Sn 
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x +t - \, Se, (42) 
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By transformaticn to Laplace netation the equations become 








k } n 
5 Saka) = = Sx(s) +2 Sn(s) 2, 2 Snls) ae Sxls) 


_-B° 
= Sn(s) + sr dy Se, (s) (43) 
ger {=1 
k n 
Se,(s) = ae S n(s) + Bite SH0 ~ A, Sey(s) (ib) 
ple~ pie” 


Solution of Equation &s for Se,(s), substitution in Fquation 


43, and rearranging yields the transfer function 
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5 6 448 
The substitution of 2 fe for 2. : in Nquation &5 
) ie, *~ 


yields 


n Le- + 


L j , 
, 7 * ayo i? ( = 735 


4 i 


The multiplication of numeratcr and deneminater by s + A and 


collection of terms yields 





on » a pase A. | oe ds 
© 52(Lyi(d + ——2 + 6)5 + 1 (——2) 
-», t hi i 


Division of the numerator and Cenominator by £*1K, yields 





k A 
sO ee ee nremaerei (48) 
k(s) ~& ~ Lan k — Lok . 
— s*+ (1+ 72 +S 2)s + AC n° 


This equation is seen te have the same denceminator as “‘quation 
_ Lak 

35 which simplifies and factors to (s +A)(s + = In the 

numerator the factor (1-4) is seen to be nearly equal to 


one when 0.0064 is substituted for /#. Fquation 48 then be- 


comes 
ons no/ 4+ 
S kts) ~ io, (49) 
s + ra 


which is the final form of the transfer function. This result 
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is the same as Equation 32 where the scurce strength was 
varied except for the censtant rultiplier tity/ L* the average 
decay constant X and fraction cf delayed neutrens drop out 
of the final expression as before, showing that the transfer 
funetion for a subcritical assembly is independent of the 
delayed neutrons and dependent on the average neutron life- 


time. 


Modulus of the Transfer Function 





l-k 
Division of the transfer functicn, Equation 49, by = 
and substitution of 5 = jw, yields 
Se Ro/i-k, (50) 
5 kts 1 + 4a 
1-k, 


The resilt cf normalizing Iquation 50 and setting it equal 


to ¥(4jw) is 


§n(s)/n, ; 
¥( jo) * SUSI/I-K, “Ty dae Lok (51) 
pe ‘oO 


It is seen from iquation $1 that the prompt neutron break 
frequency, ,, is equal te of for a subcritical assembly. 
W effectively determines k for the subcritical assembly if 
A is known, and the reverse is also true. The numerator 
and denominator are multiplied by the complex ecnjuzgate, 


Ll - 1 wd yielding 
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Y(jw) i (52) 
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the modulus of the transfer function, cr the magnitude of 


¥(jo), is 


which reduces to 
J¥(4w | = ene wt ye (53) 
f i-k, 


The square of the modulus of the transfer function lx(3@)|? 
is then 


I“(jw)|* = (94) 


1 
D+ GGer 
O 
The square of the modulus cf the transfer function is the 
basis upon which the parameters of the transfer function can 


be evaluated as will be shewn in the next secticn. 
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STOCHASTIC FVALUATION OF THL TRANSFER FUNCTION 


ihe transfer function cf an assembly can be determined 
experimentally by intreducing sinuscidal variations in the 
input and measuring the phase angle and amplitude attenuation 
with frequency of the resulting sinusecidal variation in the 
output. This is usually accomplished by sinuscidally varying 
the strensth of the neutron source or the reactivity and 
notine the resulting variation in the neutron flux. Vari- 
ations in the neutron population due to the statistical 
nature of tne fission process cause a change in reactivity 
and a fluctuation of the power level of a eritical assembly. 
these random variations in the flux level may be considered 
as a summation cf Fourier series which transform into Fourier 
integrals when the period approaches infinity (2). The 
transfer function of an assembly could then logically be ex- 
tended to represent the respense of a system tc a random as 
well as sinusoidal input. 

Consider a voltage which is a random or stochastic func- 
tion of time and limited te a small range of frequencies be- 
tween f, and f,. One can determine a value which is propor- 
tional tc the average power by squaring the function, inte- 
grating, and dividing the result by the time over which the 
integration was performed. If this "averarse power" is divided 
by Af = f- Los the frequency increment which the average 


power represents, the result is the average power or mean 
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square veltage (amplitude) ner unit frequency. This is called 
the power spectral density function. The term power spectral 
density function has been extended to apply te any function 
which ecnsists of mean square signal amplitude per unit 
frequeney of the band width investigated. 

Nuclear reactions and fission processes which give rise 
to the neutron flux in a subcritical assembly can be con- 
sidered as a random variable which has a normal or Gaussian 
distribution about some mean value of frequency of occurrence. 
The power spectral density of such a random variable is the 
averace power (mean square signal amplitude) per unit 
frequency ef the bend width measured (2). It has the unique 
property of being a constant for this type of function and 
is defined as a "white ncise”". When such a randem variable 
is passed threugh a linear system the output is a random 
variable which has been attenuated by the system character- 
istics (10). 

If a system such as the suberitical assembly is excited 
by a white ncise, which has a constant amplitude at all 
frequencies, one can analyze the output of the assembly over 
a range of frequencies to deterrine the response or transfer 
function (8). It has been shown mathematically that the 
output power density spectrum of a chain reactine system is 
proportional to the square of the modulus of the transfer 


function (11). This has also been demonstrated experimentally 
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for several nuclear reactors (4, &). 

Figure 1 shews an ifeal system for power spectral density 
measurement in which it is assumed that the electronic 
components have perfect frequency response over the range of 
the investigation ané intreduce ne neise. Under these condi- 


ticns the system equation wenld be 


Gyy(M) = Gi Co) ¥(5o)|*= KixCga)|® (55) 


i 


where GC 4) input power spectral density = eccnstant. 
6, Aa) = output pewer spectral density. 
¥(iW) = system transfer function. 


ly(4@)|* 


ft 


square of modulus of transfer function. 

Since the frequency response of the measuring equipment 
is net perfect, it alsc introduces a transfer function into 
the equation. This may be neglected ani the response of the 
measuring equipment considered independent of frequency, if 
the break frequency of the equipyent is beyond the range being 
investigated (8). The final form of the system equation was 


shown to be (8) 


Gyy(w) =A + BIY(jo)/° (56) 


where A represents the level of measuring equipment noise 
and B represents the level of the assembly noise. The sub- 


stitution of !quation 54 into lquetion 56 yields 
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Figure 1. 
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Gyy! ) + . (ay? (57) 


A sketch cf Lauation 57 is shown in Ficure 2. 
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ARBITRARY 
UNITS 


LOG POWER SPECTRAL DENSITY 


FUNCTION Gyy(w) 


LOG FREQUENCY - CYCLES /SECOND 


Figure 2. ideal power spectral density function ef an 
agsentdly 
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DESCRIPTION CF APPARATUS AND ITAXPPPIAONTAL PROCEDURE 


the graphite moderated subcritical assembly (14) utilized 
in this investigation was constructed et Iowa State University 
in 1957. Figure 3 is a photograph of the subcritical assembly 
with the end cover reseved. The blocks were constructed from 
AGK grade graphite cylinders with diameters of 7.0 in. and 
6 3/8 in., and length of 60 in. The lewer nine layers were 
made up of blecks with a 6.6 inch square cross section and 
rounded corners cf radius 3.5 inches. The top five layers 
were constructed of blocks with a 6.0 inch by 5.0 inch cross 
section with rounded corners of radius 3 3/16 inches. The 
density of the graphite was 1.56 gem. per cu. em. or 97.3 Lb. 
per cu. ft. An 8¢ inch lattice was fermed by placing uraniun 
slugs in alternate channels as shown in Figure 3. The slugs 
of vranium were 1.000 in. in diameter and &.00 in. in length 
and encased in aluminum cylinders 0.040 in. thick and with 
ends 0.20 in. thick. The outside dimensions were 1.080 in. 
in diameter and 6.40 in. in length. The uranium density was 
19.0 gm. per cn. cm. cr 1186 lb. per en. ft. 

The assenbly was covered on the top and sides with a 
0.010 in. sheet of cadmium sandwiched between a 3/6 in. sheet 
of plywood and a 1/8 in. sheet of masonite. An effective 
“plack wall" for thermal neutrons was produced by the cadmium. 
The horizontal distances between the side sheets of cadmiwn 


was 60.5 in. and 62.5 in. With the extrapolation lengths 
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Figure 3. ®ront view of the subcritical assembly 


on: 








added, the effertive size of the assembly was approximately 
62 in. (side tc side) by 64 in. (front to back). 

The suberitical assembly was mounted on a pedestal as 
shown in Figure 3. The three snaces under the assembly housed 
11.5 in. deep water tanks. The water provided protection 
from the neutron scurce which was mounted in the center tank 
uncer the pile. 

the neutron scuree consisted cr five individual 
plutenium-beryllium sources in containers 1 in. in diameter 
and 1 3/8 in. in length. The individual source strength was 


6 neutrons 


rated at cne curie emitting apyroximately 1.63 x 10 
per sec. The total scurce strength was &15 x 10° neutrons 
per sec. ‘the five sources were mounted on a pedestal in the 
eentesr tank which was filled with water. The scurces were 
located on the centeriine of the assembly directly below the 
pottom layer of graphite. Five holes were cut along the 
eenter line of the front panel of the assembly at one fcot 
vertical intervals, The holes were eut to correspond to 
vaeant channels in the lattice for insertion of the oF. 
neutron detector. The lowest hele, cne foot above the bottom 
of the assembly, will be referred te as hele one, the hole 
twe feet in height as hele two and so on. 

The eleetronic circuitry for measuring the random varia- 
tions in the neutron flux is shown in Figure 4. The neutron 


probe used was a Wood Counter Laboratory BF. detector. A 
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Nuelear Chicage Ultrasealer Model 192A was used to surply 
high voltage tc the detector and te amplify the signal from 
the detector. A Tektronix Type 112 Direct Coupled Amplifier 
with a flat frequency response from DC to 2 megacycles was 
used to further amplify the signal which was fed intc a 
Khron-Hite ultra low frequency Pand-Yrass Filter, Model 330A, 
the band width of the filter was adjustable between 4 low 
eut off frequency of 0.02 eps te a high cut off frequency of 
2000 eps. Frequency adjustment was accomplished by means of 
two dials calibrated with a logarithmic scale reading direct- 
ly in eps from 2 to 20 and multiplying switehes of 1/100, 
1/10, 1, 19 and 100. A capacitanee filter in the input 
cireuit filtered out any DC component of the signal. The 
filter comprised an extremely low noise circuit, with in- 
ternally generated hum and neise less than 100 nicrovolis. 
the filtered rancem signal was then fed into a two channel 
Brush Necorder Mark II for further amplification and record- 
ing. The recording pen sensitivity ote variable from C,01 
to 10 volts per chart line, and the chart speed was variable 
from one to 125 mm per sec. This circuitry was used to in- 
vestigate the frequency range of 0.02 ers to 20 cps. An 
alternate system was used for the frequency range of 20 eps 
to 100 eps and beyond, tc facilitate data taking and 
processing. All of the components after the Nuclear Chicago 


Ultrascalar were replaced with a Hewlett-Packard Wave analyzer 
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The power spectral 


The signal was then pas 


to 0.04 eps, 9.04 eps te 0.08 
the output of the band pass filter was 
for two 


typical 


jJensity function fer the 0.02 


signal was cheeked pericdically with an 


saturation 


sed threugh the 


ent band pass such as 


cps, 0.08 eps to 0.16 


reeorded 


minutes for 


random date for 


band pass widths are shown in Figure 6. 
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Figure 5. Electronic equipment 
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Figure 6. Typical random data for representative band-pass 
widths 
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similer te thet outlined by Murphy (12). Cne hundred 
ordinates of each random trace were measured at one secon 
intervals. Fach ordinate was squared, and then the ordinates 
and their squares were summed seperately and averaged. The 
variance, 7°, was cbtained by taking the aifference between 
the average of the squared value and the square of the average 
value. This is equivalent to the average cf the square of 


the difference from the mean. 





g-® = (y2 = F2) = (y+ #)? 


vY 


The power spectral density function G_.( ) was then obtained 
by the band ynass width in cps. This 


by dividing the variance 
resuited in a value cf mean square amplitude per unit frequen- 


cy, Or mean pover per unit frequency, for the band width in- 


vestigated. 
74° mean sau plitud 
6. (0) «= mean square amplitude 
vy Af, unit frequency 


The power spectral density for the 20 eps to 100 ens 
range was obtained by squaring the reading from the wave 
ansziyzer output meter. The meter read directly in PMS power 
per unit frequency of the band width investigated. By squar- 
ing the reading, the mean square amplitude or power per unit 
frequency of the band width investigated was obtained. 

A sample ealeculation of the type used to reduce the data 


from the Brush recorder is included in the Appendix. 
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DISCUSSION OP NPSULTS 


The results ef the data reduction are presented 
graphically in Figures 7 and & with the power speetral densi- 


ty, G(W), in arbitrery units versus freqneney in cycles 


vy 
per second. The data were plotted only from 0.03 cps te 
10.0 eps because the power spectral density appeared to be-~ 
come constant by 10.0 ens and remained a ecnstant from 10.0 
eps to 100 eps. This constant level represents the "A" term 
in Equaticn 57 which is the white noise of the instrumenta- 
tion which is due primarily tc the BF, detection chamber. 
The instrumentation noise decreases with increasing height 
in the assembly. This is explained by the fact that west of 
the instrumentation noise is due te the flux level or count- 
ing rate in the detection chamber which decreases with in- 
creasing height. 

Figure 7 shows the change in power spectral density 
function along the axis of the assembly. For hele one (z = 
1 ft) the power spectral density was determined te be ap- 
proximately a constant throughcut the range of frequencias 
investigated. This indicated that the randem variations in 
the neutron flux were primarily due te the flux from the 
source which is "white noise". Hole two (z = 2 ft) shows a 
deviation from the constant power spectral density function 
of the instrumentation and indicates an inereasing contribu- 


tion from the fission flux and external souree which have 
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Figure 8. Comparison of measured with calculated power 
spectral density function 
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been attenuated by the transfer funeticn of the assembly. 

Hole three (z = 2? ft) is 91.4 om or approximately two diffusion 
Lengths above the external source. the power spectral density 
funetion of the neutron flux at this level indicates that the 
transfer function of the assembly has cperated on the flu 

to @ much greater extent. Hole four (2 = § ft) shrews a 

maximum slope or attenuation of the power spectral density 

by the transfer fumction of the assembly. Hole five (2 = 5 

ft) shows a decreasing slepe an’ a tendency to appreach a 

break frequency teward the lower frequencies. This might be 


ez 


® Ma a mm 
: tie “2S” coaerric 


be 


explained oy 2 erease i 


.-— 


ont of Equation 


§7 which represents the level of assembly neise and is pro- 


ay 


porticnal tc the square reot of the power or flux level in 


‘3 


the assembly. As “5” decreases, the instruecntation neise 
Level approaches the assembly ncive Level, and the rounding 
off of the breeks in the eurve tends to mask the true slope 
ef the powsr spectral density function. 

He definite break frequency was diseernible from the 
Gata presented in Figure 7. ‘the data fer heles three, four, 


and five were replotted on Firsure © aleng with a plet of the 


= 


L 


theoretical sgquered modulus cf the transfer funetion frem 
Pquation 54. The theoretical equation was pletted using 
representative values of (= 9.08 see and k = 0.5. Comparison 


with the thecretical curves shows the slope of the pewer 
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spectral density of hcle four to be arprovimately the same as 


thet for the theoretiesl curve. 
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CONCLUSIONS 


the results were inconclusive because the neutron vreak 
frequency was not determinable and as a result it was not 
possible to evaluate the ratic > exactly. 

ihe best agreement between theory and experimental 
results was cbtained st zg = 4 ft = 122 em. At this position, 
the slope of the attenuated rewer spectral density function 
WAS approximately the same as that predicted by theory. The 


value of the ratio isk was determined te be approximately 


0.4/0.08 = 5.0. ni tedies indicate that this ratio changes 
as 2 increases in value to beyond two diffusion lengths where 
it remains feirly constant until masked by the white noise 
of the detection chamber. 

the evaluation of the transfer functicn of the sub- 


eritical assembly by analysis of the power spectral density 


shows promise and sheuld be pursued further. 


Mem fever «fl! eure eoieel sewers! som fv lies ast 
foe eee 18 é Greet « be Oe adisubeh ow tee cen teamed 
Wiraare 334 atyex ni eherteme ct » Olbmeey 

etme tharas Sem eared reeetes tromroane Teed of) 
-Ooldiniy sits oA om TEL aN dw 2 ie tecieids com at teen 
ontisant “leas feieergs cavcn iecemmisip eds “Ls wgole ot 
WAT ereeds wi hesntoon Fh1? ae seme on yiserabroneys mew 
COvremixetyes od of soatevetes wae Sd etter nit Ye solar 
neatady ofiet otis sade otedinar védnnd Wat Auk = 
OMY AM gre I GniaeTttt oat Seed AD ay ter a! reranvenil « 
eRe S9LMy Ady We teclnen [Jen Seetnume LyIA) seddeees 

=dgt 90) Jo miticoe? reteset? aff To wollmifeay off 
Vineah ipevesae wero acy Jo ateclene wr Udweewe Leaked 
fete feywepg oc igor) Aad wd dewey 
























7 


SUGGTI@TIONS POA PONTHOR SIUDY 


The plutorium-beryliium neviren sources utilized in this 
investigation were mounted on a pedest2i beneath the assembly 
and were unmoderated except for the moderation by the graphite 
of the assembly. Interesting results wight be cbhtained by 
enceasing tne sources in a bloek of paraffin moderator, and 
surrounding the moderated scurces with graphite instead of 
the water shield. Nore signifiernt results might be obtained 
by the use of a mechanical or electronic device for squaring 
and averacing the data. This means of data processing would 
allow much more rapid and accurate determinstion of the power 
spectral censity, hence breadening the scope ef investigation 


possible. 
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